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ABSTRACT 

We present the rest-frame optical galaxy merger fraction between 0.2 < z < 1.2, as a function of 
stellar mass and optical luminosity, as observed by the Canada-France-Hawaii Telescope Legacy Deep 
Survey (CFHTLS-Deep). We developed a new classification scheme to identify major galaxy-galaxy 
mergers based on the presence of tidal tails and bridges. These morphological features are signposts 
of recent and ongoing merger activity. Through the visual classification of all galaxies, down to 
ivega < 22.2 (sa27,000 galaxies) over 2 square degrees, we have compiled the CFHTLS Deep Catalog 
of Interacting Galaxies, with « 1600 merging galaxies. We find the merger fraction to be 4.3% ±0.3% 
at z ~ 0.3 and 19.0% ± 2.5% at z ~ 1, implying evolution of the merger fraction going as (1 -|- z)"^, 
with m = 2.25 ± 0.24. This result is inconsistent with a mild or non-evolving (m < 1.5) scenario at a 
>4cr level of confidence. A mild trend, where massive galaxies with Af* > lO^^'^M© are undergoing 
fewer mergers than less massive systems (M» ~ IO^^Mq), consistent with the expectations of galaxy 
assembly downsizing is observed. Our results also show that interacting galaxies have on average 
SFRs double that found in non-interacting field galaxies. We conclude that (1) the optical galaxy 
merger fraction does evolve with redshift, (2) the merger fraction depends mildly on stellar mass, with 
lower mass galaxies having higher merger fractions at z < 1, and (3) star formation is triggered at all 
phases of a merger, with larger enhancements at later stages, consistent with A^-body simulations. 
Subject headings: galaxies: evolution - galaxies: formation - galaxies:interactions - galaxies: starburst 



1. INTRODUCTION 

In the past three decades the concept of galaxies as "is- 
land universes" slowly evolving in isolation has changed 
dramatically. Gravitational interaction between galaxies 
is now considered a relevant factor in a galaxy's evo- 
lution, capable of altering its morphology, luminosity, 
color, size, star formation rate (SFR), and mass distri- 
bution, all over a relatively short timescale. According 
to the A-dominated cold dark matter (ACDM) model, 
the merger rate of dark matter halos and similarly the 
merger rate of galaxies are the most fundamental pro- 
cesses in structure formation. Numerous A'^-body simula- 
tions and semi-analytical models have studied the merger 
rate of dark matter halos, predicting ev olution with red- 
shift as (1 + z)", with 1.0 < m < 3.5 (iGottlober et al.l 
[200lHBerrier et al.|[200llFakhouri fc Mall2008[ ). A direct 
comparison between dark matter halo merger rates and 
galaxy merger rates is difficult due to the uncertainty in 
the galaxy halo occupation number. Although measuring 
the frequency that galaxies merge as a function of cos- 
mic epoch can place powerful constraints on the theory 
of galaxy evolution and structure formation. 

Traditional observational approaches aimed at investi- 
gating the major galaxy merger rate, measure the fre- 
quency of galaxy mergers or galaxies in close pairs, 
spanning a range of redshifts. The merger or close 
pair fraction should evolve in a similar manner as the 
merger rate assuming the timescale over which the 
merger selection criteria are sensitive to is indepen- 



dent of redshift. iToomrd (|1977[ ) was the first to sug- 
gest that the merger rate may be larger at higher red- 
shifts, by using estimates of past merger remnants. 
Numerous studies have estimated the evolution of the 
merger rate, yielding highly discrepant measurements 
of m, ranging from no evolution (m ^ 0) to strong 
evolution (m ^ 5) (iZepf fc Kool 119891: iCarlberg et al.l 
1991 [2OOOI: iPattonet al.l Il997t iLe Fevre et al.l | 200C : 
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Cpnsehce et al."2003' 'Lin et al."2004'; 'Bri dge et aL|[2007 
Kartaltcpc et al. 2007; Lotz et al. 2008a, to name a few) 
. It is often suggested that this discrepancy stems from 
the variety of techniques used to identify galaxy merg- 
ers, incompleteness corrections, as well as differences in 
sample selection. Although these factors likely play a 
role, typical sample sizes have been on the order of a 
few ten's of mergers (excluding Kartaltepe et all ()2007D ) 
resulting in > 30% uncertainty in the evolution of the 
merger rate and the effect of cosmic variance has been 
largely unexplored. 

Over the past decade it has become clear that the 
average SFR per unit comoving volume (CSFR) has 
declined by an order of mag nitude since z ~ 1 
(jLillv et all 119961: IMadau et alTll998D . A fundamen- 
tal question that remains is the cause of the drop of 
the CSFR. It has been suggested by some close pair 
(iBurkev et al."1994l : lYee fc Elling son"199 5l: iPatton et al.l 
|1997, 2002; Bridge eTaLl 120071: iK artalt epe et al.l |2007D_ 
and merger studies ^ Lc Fevr e et al.li2000tlConselice et al.l 
120031 : IKampczvk et al. ,2003) that an increased merger 
rate at higher redshift is either partially or completely 
responsible for the higher SFR density at z ~ 1. How- 
ever, there are merger studies wh ich have found little 
or no evolution of the merger rate fCarlberg et al.l j2000l : 
IBundv et a l. 2004; Lin ct al. 2004; Lotz et al. 2008^), in- 
dicating the decrease in the volume-averaged SFR den- 
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sity since ~ 1 is not a result of a declining merger 
rat e but rather decli ning star formation in disk galax- 
ies (|Lotz et al.ll2008aD . 

In order to probe the merger rate a clear definition 
of what constitutes a merger and a robust identifica- 
tion technique are required. One traditional method is 
to search for close galaxy pairs. However, even with 
radial velocity measurements for both galaxies, about 
half of all physical pairs may be chance superpositions 
(jPatton et al.il2000l I200I . Another method is to calcu- 
late internal asymmetries, but these quantitative mor- 
phological parameters are subject to complication, since 
star formation itself is sufBciently violent and chaotic 
and therefore not smoothly distributed, potentially mim- 
icking features of a merger when none has occurred. 
However, long tidal tails are nearly a foolproof signa- 
ture of imminent mergers. Moreover, tidal tails are 
a relatively simple, completely dynamical phenomenon 
that can be studied in de t ail with A^-body simulations 
(iToomre fc Toomrd 119721: iBarnesI 119921: iHopkins et all 
|2008() . In spite of these positive features little work on 
galaxies with tidal tails has been done at high redshift. 
The reason is simple: tails have a relatively low surface 
brightness. 

The Deep component of the Canada France Hawaii 
Legacy Survey (hereafter CFHTLS-Deep) provides some 
unique advantages for studying morphologies and the 
merger rate of galaxies. The deep ground-based imag- 
ing {i' < 27.3) is ideal for the detection of low surface 
brightness features like tidal tails. The survey area is 
also spread over four fields allowing for cosmic variance 
estimates of the merger rate, while the larger area al- 
lows for merger fraction estimates no longer dominated 
by small number statistics as seen in previous studies 
(excluding Kartaltepe et al. 2007). 

In this paper, we apply a new technique for identifying 
major merger candidates in the CFHTLS-Deep Survey. 
Through the visual classification of «27,000 galaxies over 
2 square degrees, we have complied the CFHTLS-Deep 
Catalog of Interacting Galaxies^ which contains wlGOO 
galaxies between 0.1 < z < 1.2 with tidal tails and 
bridges. This is the largest catalog of its kind in the liter- 
ature. In Section 2, we describe the optical photometry 
from the CFHTLS-Deep Survey, along with the method- 
ology used for deriving the photometric redshifts, stellar 
masses and SFRs. Section 3 outlines the technique used 
to identify interacting galaxies, including completeness 
tests. We measure the galaxy interaction fraction in Sec- 
tion 4, its dependence on optical luminosity, stellar mass, 
and address cosmic variance. With additional assump- 
tions about the merger timescale we estimate the merger 
rate and interaction history of galaxies from z ~ 1 in Sec- 
tion 5, followed by the impact mergers have on the SFRs 
of galaxies (Section 6). We summarize our conclusions 
in Section 7. 

In the discussion that follows, we assume Vega magni- 
tudes and any calculation requiring cosmology assumes 
r^M 0.3, r^A = 0.70, and Ho=70kms-i Mpc^i. 

2. THE CFHTLS DEEP OBSERVATIONS 

2.1. MegaCam Optical Imaging 

The CFHT Legacy Survey (CFHTLS) is a joint com- 
munity project between Canada and France with more 



TABLE 1 
CFHTLS-Deep: Final Stacks 



Field 


Filter 


Int. Time [hrs] 


Limiting Magnitude 




Dl 


u* 


10. 6h 


26.7 


(0.96 square degree) 


g' 


9.5h 


27.1 




r' 


18.8h 


26.8 




i' 


45.1h 


26.4 




z' 


20.0h 


25.6 


D2 


u* 


3.4h 


26.4 


(0.90 square degree) 


g' 


5.7h 


27.0 




r' 


10.7h 


26.5 




i' 


22.2h 


26.2 




z' 


12.0h 


25.2 



Note. — Limiting magnitudes are estimated by adding artifi- 
cial galaxies, faded and rescaled into the images and then trying to 
recover them. 



than 450 nights over a 5 year period that commenced in 
2003 June. The survey exploits the square degree Mega- 
Cam camera (Bouladc et. al. 2003) on CFHT which has 
36, 2048x4612 pixel CCDs with a pixel scale of 0".187. 
The Deep survey is one component of the CFHTLS 
(deep, wide, very wide), covering four low Galactic 
extinction fields in u* , g' ,r' ,i' , z' . Each survey field, 
named Dl through D4 is one square degree and dis- 
tributed in right ascension (R.A.) for efficient observ- 
ing throughout the year. The survey fields were cen- 
tered at J2000 R.A.= 02^26'"00^ decl.= -04'i°g30'00" 
(Dl), R.A.= 10''00™29^ decl.= -K02'i'=gl2'2l" (D2), 
R.A.= 14''17"54", decl.= +52'i°s30'3l" (D3), R.A.= 

22'*15"31", decl.= -17'^'=S44'05" (D4). The analysis 
that follows utilizes only the Dl and D2 fields, a total 
of 2 square degrees, due t o their extensive ancillary data 
(see lSuUivan et al1l2006al for details). 

The optical images used to derive the galaxy param- 
eters, and morphological classification we re constructed 
by the Supernova Legacy Survey (SNLS; iSullivan et al.l 
2006b), a key project in the CFHTLS. Deep optical 
stacks were generated for each filter {u*,g', r' , i', z') with 
an imposed wavelength-dependent seeing limit (~4 pix- 
els) to ensure a high degree of resolution. The goal was 
to maximize exposure depth while retaining similar and 
excellent seeing in the different filters. The typical seeing 
of the final stacks is 0".7-0".8 [i' band). The individual 
"Elixir" processed images are available from the Cana- 
dian Astronomy Data Centre (CADC); however, the fi- 
nal stacks and photometric catalogs including redshifts 
in which this work is based are not currently public. 
The details of the processing technique are discussed in 
([Sullivan ct al. 2006b). 

Source extraction and photometr y were performed 
on ea ch Deep field using SExtractor (jBertin fc Arnouti 
IT996h in dual image mode. Detections were performed 
in the i' filter («' ^ 26.3) and photometry measurements 
calculated in each of the five filters, u* , g' ,r' ,i' , z' . The 
exposure times and limiting magnitudes reached for each 
filter in both fields are given in Table [1] We applied a 
bright star and bad pixel mask to the images prior to 
running the source detection to eliminate noisy or con- 
taminated regions. The total area masked is < 10% for 
each field. We compared our number counts to those 
from the Hubble Deep Field (HDF) North and South 
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(IWillia.ms et alJ [l996l: iMetcalfe et alJ [200l and COS- 
MQS (jLeauthaud et al.ll2007D . Thev were found to be 
consistent for galaxies brighter than i' < 24, on which 
this work is focused. 

2.2. Galaxy Properties: Redshift, Stellar Mass, SFRs 

In this section, we outline the methods used to con- 
vert the optical fluxes of our sources into a photometric 
redshift, and derive other properties such as stellar mass 
and SFR. 

In order to study the potential evolution of the in- 
teraction and star formation rates of galaxies, we need 
to derive a redshift estimate for each galaxy. Although 
spectroscopic redshifts are the most precise distance mea- 
sures, they are observationally expensive to obtain for 
large samples. A strength of the CFHTLS-Deep survey 
is its high quality five-band optical imaging, which can be 
used to derive a photometric redshift estimate. The five 
optical bands can be combined to produce a broadband 
spectral energy distribution (SED) that can be compared 
to a set of template SEDs to estimate properties of the 
:alaxy such as its redshift, age, stellar mass, and SFR 
Baumlll962t iLoh fc Spillail 119861 ). This technique uses 
broad spectral features like the 4000A discontinuity or 
the Lyman break for comparison to templates. 

We fit a series of galaxy template SEDs to the broad- 
band fiuxes of each galaxy. The best-fit SED is deter- 
mined through a standard minimization procedure 
between the synthetic photometry generated by inte- 
grating the template SEDs through the CFHTLS fil- 
ters, and the observed fluxes (including the flux errors). 
This was perforrned with the Z-Peg template fitting code 
(iBolzonella et al.l 120001 : iLe Borgne fc Rocca-Volmerangd 
Its methodology is similar t o that applied by 
other photometric redshift codes fe.g. lGwyn fc Hartwickl 
[19961) . 

We employ a set of synthetic templates 
computed with the P EGASE-ii galaxy evo- 
lution code (iFioc fc Rocc a-Volnicran g3 119971: 



Le Borgne fc Rocca-Volmerangei l2002t ILe Borgne et al.l 
20041 ). Both the SEDs and the photometric redshift 



code have b een extensively tested and used in the lit- 
erature I'e.g. iMcCarthv et al.|[20p: iGrazian et al.l 120061: 
ISullivan et al.ll2006bD . We use eight evolutionary sce- 
narios that evolve self-consistently with age and assume 
a Rana-Basu initial mass function (Rana fc Basu 1992). 
These scenarios were designed to match the average 
colors of local galaxies and to reproduce deep galaxy 
number counts. Considering several tens of time steps 
for each scenario, the total number of synthetic SEDs 
is «500. When fitting a given galaxy only templates 
younger than the age of the universe at the redshift of 
the galaxy are considered. The accuracy of the photo- 
metric redshifts is determined by comparin g them to the 
SNLS spectroscopic samp le in our fields ([Howell et al.l 
120051: iBronder et all 120081) . The photometric redshift 
accuracy down to i' ^ 22.5 is (TAz/i^ + Zg) = 0.04. 

Two additional physical parameters of particular inter- 
est are the SFRs of interacting galaxies and their stellar 
masses. T hese quantities are both derived using the Z- 
Peg code (ILe Borgne fc Rocca-Volmerangd [213021) . The 
mean recent SFR for a galaxy is determined from the 
best-fit SED, averaging the SFR over a period of 0.5 



Gyr. The total stellar mass of a galaxy was derived 
by integrating the total star formation history (SFH) 
of the best-fit scenario, up to the best-fit age and sub- 
tracting off the mass from stars that have died. The 
results from these techniques were found to be in good 
agreement with template fits when a spectroscopic red- 
shift was used ([Sullivan et al.l [2006bl) . The systematic 
errors associated with these techniques, their application 
to CFHTLS d a ta, and consistency checks are outlined in 
ISullivan et aP ([2006bD . 

3. CFHTLS-DEEP CATALOG OF INTERACTING 
GALAXIES 

3.1. Identifying Interacting Galaxies 

The first step in studying the frequency of galaxies 
involved in an interaction is to define a clean, robust, 
and useful definition of an interacting galaxy. Morpho- 
logically, interacting galaxies can exhibit long tidal tails, 
bridges (linking two or more galaxies), ring structures, 
stellar bars and/or enhanced spiral structure and fre- 
quently appear severely distorted. We have focused our 
identification methods on confidently selecting galaxies 
which have recently undergone or are presently under- 
going a tidal interaction. The presence of a tidal tail or 
bridge is incontrovertible evidence of a recent interaction. 

There are two different avenues one can take to mor- 
phologically select interacting galaxies. The first is to uti- 
lize quantitative morphological software that measures a 
galaxy's structural param eters, su ch as asymmetry, Gini 
or M^n ([Abraham et aLi .1996b, al: iConselice etall [2OOOI : 
iLotz et al.ll2004[ ). As a galaxv undergoes a merger, tidal 
fields distort the galaxies radially, drawing out galactic 
material into long plumes and tails. These structural 
parameters can help describe the level of a galaxy's dis- 
turbed appearance and are commonly used to identify 
mergers. 

This is an efficient and automated approach. How- 
ever, a complication of this method is the need for high- 
resolution data, especially at high redshifts when Hubble 
Space Telescope (HST) observations are a requirement. 
A further complication is that HST /ACS images do not 
easily detect low surface brightness features, like tidal 
tails due in combination to the small plate scale and point 
spread function (PSF). 

The second method of identifying interacting galaxies 
is based on pure visual inspection. This technique has 
been utilized by many in the past to morphologically clas- 
sify galax ies (Hubble 1926; Vorontsov-Velvaminovl [19591 : 
lSandagd[r961- ,Arp,,196& to name a few). Qualitative 
classification is able to identify the low surface bright- 
ness features that the automated method has difficulties 
with, as it utilizes one of the best pattern recognition 
computers-the brain-eye combination. A drawback of 
this method is that it can be laborious, subjective and 
can suffer from reproducibility issues when the person 
conducting the classification changes. However, when a 
set of visual criteria are clearly defined, and the features 
themselves are striking (e.g., long tidal tails and bridges), 
visual inspection can be a highly accurate method of mor- 
phological classification (see Section [3T3| . 

The ground-based nature of our dataset coupled with 
its depth make it ideal to detect low surface brightness 
features down to <29 mags arcsec"^. Ultimately, we 
define an interacting galaxy to be one with a tidal tail 
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or bridge. These merger signatures are visible after a 
first encounter and can persist even after the galaxies' 
nuclei have coalesced. Although our primary objective 
was to identify galaxies undergoing an interaction, other 
classifications were included and will be presented in a 
future paper. 

The morphology of galaxy interactions relies on a large 
number of variables, such as the geome try of the en - 
counter, impact velocity, and mass ratios (IBarneslll992l) . 
By defining a sample of tidal tailed galaxies we are se- 
lecting galaxy mergers which have recently experienced 
a merger of mass ratio M1/M2 > 1/10. Galaxy en- 
counters with mass ratios > 1/10 have been s hown to 
have significant impact on a ga laxy's evolution ()Barnesl 
[T99I IMihos fc HernauistI 119961) . Although some of the 
tailed galaxies in our sample will have been a result 
of lower mass ratio or minor mergers (defined as those 
with 1/10 < M1/M2 < 1/4), the majority are likely 
to have been the consequence of M1/M2 > 1/4: inter- 
actions, due to the marked nature of the tidal features 
before being includ ed in our sample (e.g. [Barnes! [19921 : 
IMihos fc Hernauistlll996l:IDubinski et al.lll999D ^ 

3.2. Classification Methodology and Scheme 

The large area and deep optical imaging of the 
CFHTLS-Deep survey makes it well suited to identify- 
ing interacting galaxies. All galaxies in the Dl and D2 
fields (2 square degrees) down to ivega = 22.2 (w 27,000) 
were visually classified using the survey's deepest filter, 
the i' band. The apparent magnitude distributions were 
similar for the Dl and D2 field. There were small dif- 
ferences (< 10%) at the bright end but were due to the 
Dl field covering 0.05 square degrees more than D2, cou- 
pled with a larger number of bright stars in the D2 field, 
resulting in a larger masked out region. 

Eac h galaxy was yiewed on a computer screen with 
DS9 flJove fc Mandell [2003D . and the classification was 
logged through keypad entry to a file. More than 90% of 
galaxies with apparent magnitudes brighter than i' <21.9 
were able to be visually classified. Although some galax- 
ies at these apparent magnitudes had angular sizes too 
small for a reliable classification to be made, the ma- 
jority of the remaining 10% lied within halos of bright 
stars, or bad regions of the images that were not masked 
out. At faint magnitudes visual classification becomes 
increasingly difficult and less reliable, due typically to 
the smaller angular size of the objects. We therefore im- 
posed an apparent magnitude cutoff of i'y^ga <21.9 to 
ensure highly reliable classifications. Our final morpho- 
logical catalog consists of 25,194 galaxies. 

Galaxies were classified as interacting or merging 
if they exhibited a tidal tail or tidal bridge between 
galaxies. These two tidal signatures encompass early 
through to advanced stage mergers. We then further 
sub-classified these objects in a simplistic way according 
to the length of their longest tidal tail as well as the 
number of tails to see if they were correlated with 
other properties such as the SFR or active galactic 
nucleus (AGN) activity. The length of the tidal tail was 
defined in relation to the angular size of the host galaxy. 
Figure 1 provides i' band and color composite {g' ,r' 
examples of galaxies classified as interacting. Galaxy 
mergers were identified and classified into the following 
primary categories: 



Long tidal tailed galaxies - Galaxies exhibiting a tidal 
tail longer than the diameter of the host galaxy. These 
mergers primarily represent intermediate to late stage 
interactions after a first encounter. 

Medium tidal tailed galaxies - These systems have a 
tidal tail length approximately equal to the diameter of 
the host galaxy, and also probe interactions after the 
first passage to advanced mergers. 

Short tidal tailed galaxies.- This classification iden- 
tifies galaxies with tidal tails that are less than the 
diameter of the host galaxy. This classification suffers 
from the most contamination from spiral arms being 
misidentified as tidal tails. Therefore, short tidal tailed 
galaxies are only included in the potentially interacting 
sample and not considered in the primary analysis of 
this paper. 

Tidally bridged galaxies- Galaxies found to have 
a tidal bridge linking it to another galaxy. This 
classification represents a narrower time frame in the 
merger sequence, typically immediately after the first 
encounter, or second passage, and is in turn well suited 
to comparative studies with simulations. 

Double nuclei.- Our final primary merger classification 
identifies mainly advanced stage mergers, with two or 
more nuclei overlapping in a common envelope. Since 
some galaxies with apparent double nuclei are a result 
of line-of-sight projections this classification was further 
separated into those galaxies with a multiple nuclei and 
additional evidence of a tidal tail and those with no 
tidal tail signatures. 

Table [5] outlines the classification scheme and the num- 
ber of galaxies identified in each merger class. Ulti- 
mately, we find 1586 galaxies with clear signatures of 
a tidal tail or bridge, with an additional 970 potentially 
interacting galaxies. 

The sample of potentially interacting galaxies is com- 
prised of galaxies with a double nuclei with no evidence 
of a tidal tail, which can be a result of projection, and 
those with "short" tidal tails since these tails are less pro- 
nounced and can be confused with spiral arms. However, 
if the "short" tailed galaxy had secondary evidence, such 
as the presence of a double nuclei or a tidal bridge they 
were considered to be undergoing a merger. Since an aim 
of this paper was to measure a lower limit of the frac- 
tion of galaxies undergoing a major merger the sample 
of potentially interacting galaxies was included only as 
secondary measure of the merger fraction (see Section [4]) 
and not in the merger rate determination or the general 
conclusions of this work. 

3.3. A Re- Classification Experiment 

A potential bias in using a visual classification tech- 
nique is its reproducibility when other individuals in- 
spect the same sample of galaxies. A key criteria of this 
work, however, was the requirement of strong tidal sig- 
natures before a galaxy was deemed "interacting" , dra- 
matically reducing the classification variance by different 
individuals. C.B. visually classified all galaxies under 
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Fig. 1. — Examples of galaxies classified to have tidal tails or bridges from the CFHTLS-Deep Catalog of Interacting Calaxies. The 
i' band images (left) and color composites {g',r',i') range in class from galaxies with tidal tails (various lengths), close pairs with tidal 
bridges to double nuclei with tidal tails. Each stamp is 100/i~^kpc on a side. The white circle marks the galaxy that has been classified. 
The XY a,xes are in arcseconds. 
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TABLE 2 

Classification Scheme for Interacting Galaxies 



General Description 


JV 


"Long Tidal Tails" : Tail Length > Diameter of Host 


Galaxies with a long tidal tail: 1 tail 
Galaxies with a long tidal tail: 2 tails 
Galaxies with a long tidal tail: > 2 tails 


157 
183 
36 


"Medium Tidal Tails" : Tail Length ~ Diameter oi Host 


Galaxies with a medium length tidal tail; 1 tail 
Galaxies with a medium length tidal tail; 2 tails 
Galaxies with a medium length tidal tail; > 2 tails 


293 
138 
25 


Close Galaxy-Galaxy Pairs with Tidal Bridges 


Galaxy in a close pair with bridge, no tail 
Galaxy in a close pair with bridge, + short tidal tail 
Close Pair with bridge + medium tail(s) 
Close Pair with bridge + long tail(s) 


486 
48 
49 
50 


Double Nuclei (DN) 


Galaxies w/ a DN + short tail(s) 
Galaxies w/ a DN + medium tail(s) 
Galaxies w/ a DN + long tail(s) 


12 
47 
62 


Total Number of Confidently Interacting Galaxies 


1586 


Potential Galaxy Interactions 


"Short 'I'idal I'ails" : Tail Length < Diameter oi Host 
Galaxies with a double nuclei (DN), no tidal tail 


586 
384 



Note. — A'^ specifies the number of galaxies in each morphological 
classification 



study. To cross-check these classifications, M.S. and K. 
Bundy (KB) classified a sub-sample of 700 objects. This 
sample of galaxies was randomly selected from within 
four weighted parent classifications (490 galaxies with 
tidal tails, 110 galaxies with a tidal bridge, 50 galaxies 
classified as potentially interacting, and 50 deemed non- 
mergers). 

This re-classification experiment was a blind study, 
both K.B. and M.S. were given no information regarding 
the galaxies they were classifying. The classification was 
performed in the same manner as described in Section 
13.21 Of the galaxies used in the classification experi- 
ment, 86% had been categorized as "interacting" prior 
to the experiment (by C.B.). Both K.B. and M.S. clas- 
sified 88% and 87% of the galaxies in the experiment to 
be "interacting". The additional l%-2% found by K.B. 
and M.S. to be interacting were predominantly classified 
as "potential" galaxy interactions with < 0.5% deemed 
non-mergers by C.B. The strong agreement between dif- 
ferent individuals regarding which galaxies are "interact- 
ing" is a result of the robust visual criteria a galaxy must 
exhibit before it is classed as undergoing an interaction. 
We are therefore excluding many potentially interacting 
systems as well as those with tidal features below the 
sensitivities of our survey. The results that follow are 
therefore secure lower limits. 

As an additional test, C.B. reclassified a set of 500 
galaxies (350 interacting, 150 non-interacting) randomly 
selected. This blind self test addresses the reproducibil- 
ity of the authors own classifications. The classifications 
remained the same 97% of the time. The 3% variation 
was primarily a result of galaxies being classed into mor- 
phologically similar types, for example, an intermediate 
tidal tailed galaxy being classed as having a short tidal 



tail, and not into grossly different types. 

3.4. Classification Completeness 

When investigating how the fraction of interacting 
galaxies changes from low-z to high-z, it is important to 
investigate how tidal features become less resolved and 
fainter due to cosmological effects. To address this ques- 
tion we have artificially redshifted bright nearby galaxies 
with tidal features in our sample out to higher redshifts. 
We then reclassify the redshifted galaxies to determine 
the point in redshift space that the tidal signatures are 
too weak, and the galaxy would not be classified as in- 
teracting. 

When simulating an image of a low redshift galaxy 
to how it would appear at a higher redshift, multiple 
factors must be considered. First, we rebin the image, 
ensuring that flux is conserved, to account for the re- 
duction in apparent size of the galaxy when it is viewed 
at higher redshift. Second, the dependence of surface 
brightness (SB) with redshift (SBoc (1 -1- z)~^), as well 
as fc-correction effects needs to be accounted for. In order 
to accurately simulate the appearances of galaxies at high 
redshift we carefully selected the lowest redshift galaxies 
possible with similar luminosities to those probed at the 
high redshift end of our data {z ~ 0.7 — 1.0). Our resul- 
tant sample consisted of 54 galaxies with tidal features 
between 0.3 < z < 0.45 and Mg < -21.0. 

At redshifts between 0.4 < z < 0.9, the i'-band, which 
was the filter used for visual classification, is probing 
rest frame g'. Therefore, to help account for the fc- 
correction the r'-band images of our sample of low-z 
galaxies (0.3 < z < 0.45) were artificially redshifted, 
as they also probe rest g', like that of the i'-band im- 
ages used for classification of the main sample. When 
redshifting the low- 2: galaxies to our highest redshift bin 
(z ^^1-1.15), ideally the g'-band images should be used 
so that the same rest-wavelengths are being compared in 
original and redshifted samples. However, the g' images 
were not sufficiently deep, and in turn the r'-band images 
were used. 

Ultimately, a sample of 54 low redshift galaxies were 
redshifted out to five different redshifts, z =0.55, 0.70, 
0.85, 1.0 and 1.15. The galaxies were re-classified at 
each interval for the presence of tidal features. The re- 
covery rate remains high out to z ~ 0.85, where we can 
still identify strong tidal features in ~ 80% of the red- 
shifted galaxies. Throughout the analysis, the recovery 
rate (75% a.t z ~ 0.9) is referred to as a completeness 
limit, as it quantifies our survey's sensitivity to tidal fea- 
tures. The rate of recovery could also be used as a cor- 
rection factor when determining the fraction of galaxy 
mergers at a given redshift (as discussed in Section 2]) . 

3.5. Sample Selection 

In order to trace the fraction of merging galaxies with 
redshift, we first need to define a sample of galaxies that 
probe the same stellar mass or luminosity ranges over 
the entire redshift range. In Sections 13.21 and 13.41 we 
found that visual classifications could be securely made 
to ivega < 21.9, with > 75% completeness out to z = 0.9. 

The majority of merger rate studies to date impose 
luminosity rather than stellar mass limits on thei r data 
sets (excluding iConselice et al..,2003, ,2008; .Bundv et al.l 
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120091) . A potential bias with this selection criteria is the 
inclusion of interacting galaxies into the sample which 
prior to the merger event may have been fainter than 
the luminosity cutoff. This selection bias may lead to 
higher me rger fractions o r close pair c ounts (Bund v et al.l 
[2004: Berrier et al.l [20061 iPatton fc At ficld 2008). How- 
ever, stellar mass estima .tes are less affe cted by merger 
induced star formation ([Conselicd I2006D . and therefore 
a less biased selection criteria. In the sections to follow 
we explore the galaxy interaction fraction using a stellar 
mass limit of logM^XM^) > 9.5 (see Figure[2]), however 
also consider luminosities limits to fairly compare our 
findings with previous works. The stellar mass and lu- 
minosity limits were chosen by balancing incompleteness 
with sample size. 

To first order, our sample was divided into "interact- 
ing" and "non-interacting" galaxies. As outlined in Sec- 
tion 13.11 we define "interacting" or "merging" galaxies 
as those with tidal tail features, more specifically, with 
the following visual classifications: medium to long tidal 
tails, galaxies in a close pair linked with a tidal bridge, 
and those with a clear double nuclei and a tidal tail. We 
find a total of 1586 interacting galaxies within the Dl 
and D2 fields combined. 

In order to investigate the properties of interacting 
and merging galaxies, we also need to establish a com- 
parison sample of non-interacting galaxies. We combine 
all the "non-interacting" classifications, which are pri- 
marily composed of spiral or disk dominated sources to 
construct a fair comparison sample. We have identified 
22,268 non-interacting galaxies. 

We applied a traditional template fitting technique as 
described in Section [2?2] to acquire photometric redshift 
estimates for each of our sources. The redshift distribu- 
tion for both the interacting and non-interacting popu- 
lations (Figure [3]) is found to be comparable. This pa- 
per probes the merger fraction over the redshift range 
0.1 < z < 1.2. The upper redshift limit of 1.2 is due 
in combination to the lack of sources at redshifts > 1.2 
(caused by the apparent magnitude limit required for ac- 
curate visual classification), and the fact that the 4000 A 
break begins to move beyond our bluest filter, reducing 
the accuracy of the photometric redshifts. Ultimately, 
the CFHTLS Deep Catalog of Interacting galaxies con- 
tains 1240 interacting galaxies between 0.1 < z < 1.2, 
ivega < 21.9, with stcUar masses > 

4. THE GALAXY INTERACTION FRACTION (GIF) AT 
0.1 < Z < 1.2 

Tidal tails and bridges are a result of gravita- 
tional encounters betwe en two or more galaxies (e.g. 
iToomre fc Toomri 119721 ). and are nearly foolproof sig- 
natures of a recent or ongoing galaxy interaction. The 
statistics of galaxies with extended tidal features (tidal 
bridges and tails) is a powerful tool to study the evo- 
lution of the galaxy interaction fraction (GIF). In this 
section we present our analysis of the merger fraction 
between 0.1 < z < 1.2 using the CFHTLS Deep Cata- 
log of Interacting Galaxies. This catalog is the largest 
sample of interacting galaxies in the literature over this 
redshift range. 

The number of galaxies exhibiting strong tidally in- 
duced structures like tails and bridges is a relatively sim- 
ple and robust, lower limit measure of the interaction 
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Fig. 2. — Stellar mass as a function of redshift for the CFHTLS 
Dl and D2 fields. Galaxies classified as non-interacting arc shown 
as open circles (black) and those undergoing an interaction, filled 
circles (red). The dashed line defines the stellar mass limit (log 
mass > 9.5M0) imposed. The histogram on the right notes the 
stellar mass distribution of the field. 
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Fig. 3. — Redshift distribution for all galaxies classified as non- 
interacting (upper) and interacting (lower, dashed). The mean 
photometric redshifts for the two samples are 0.51 (non-interacting) 
and 0.58 (interacting). The two samples have comparable redshift 
distributions. 

fraction. The CFHTLS-Deep survey is sensitive to tidal 
features with surface brightnesses down to i' ^ 29 mag 
arcsec~^. Our approach is simply to compare the num- 
ber of galaxies with tidal features (Nint) to the total 
number of galaxies within the same absolute magnitude 
or stellar mass ranges [Nxotai), as a function of redshift 
(see Equation[T]). In the case where a pair of galaxies are 
connected by a tidal bridge, each galaxy will be counted 
as interacting (Nint = 2), while a galaxy with a tidal tail 
and no bridge linking it to another object is counted as 
a single interaction {Njnt = !)■ 



Galaxy Interaction Fraction{GI F) 



N,,, 



Total 



(1) 



In Figure 2) we plot the fraction of galaxies with tidal 
tails, as a function of redshift, assuming i' < 21.9, and 
stellar masses > 10^'^Mq. As previously described in 
Section [3?2l the interaction classes included in this mea- 
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TABLE 3 

Inferred Galaxy Interaction Statistics 



Redshift 




^Total 


Interaction Fraction [%] 


0.10-0.25 


108 


1745 


6.2±0.6 


0.25-0.40 


160 


3725 


4.3±0.3 


0.40-0.55 


209 


4736 


4.4±0.3 


0.55-0.70 


280 


4700 


6.0±0.4 


0.70-0.85 


258 


3413 


7.6±0.5 


0.85-1.00 


168 


1497 


11.2±0.9 


1.00-1.15 


57 


300 


19.0±2.5 


Note. — 




specifies 


the number of interacting 



galaxies. Nj-otal is tlie overall number of galaxies in 
the sample. Errors are Poisson statistics. 

surement are "medium" and "long" tidal tailed galaxies, 
those in a close galaxy pair with a tidal bridge, and those 
galaxies with a double nuclei and addition morphological 
evidence of an interaction such as a tidal or bridge. 

We find that the fraction of galaxies involved in a 
merger rises with redshift. Meaning, more galaxies were 
undergoing a tidal interaction when the universe was 
about half its current age. At low redshift (2 ^ 0.3) 
the GIF was found to be ~ 4%±0.3, and nearly triples 
by z -0.95 to 11% ± 0.9 (see Table E]). We character- 
ize the evolution of the galaxy interaction fraction by 
fitting a simple power-law increase with redshift of the 
form GIF=GIFo(l -f- z)™, where GIF,, is the present 
day interaction fraction, and m is the power-law in- 
dex. When all redshift bins in our sample are included 
(0.1 <z< 1.15) we find a best fit m of 2.25 ± 0.24 and 
GIFo of 2.15%±0.25, weighted by the GIF error bars in 
each redshift bin. 

As discussed in the following section there may be some 
potential biases at the low and high redshift ends of our 
analysis. Therefore, the evolution of the GIF was also 
fit using various redshift ranges as shown in Figure 21 If 
the lowest redshift bin (0.1 — 0.25) is excluded, a best- 
fit power-law index of m = 2.80 ± 0.25 is derived. If 
the highest redshift bin 1.0 — 1.15 (the least confident) is 
excluded m = 1.95 ± 0.25, the lowest degree of evolution 
found in our analysis. Finally, if both the lowest and 
highest redshift bins are removed, a value for m of 2.56 ± 
0.24 is found. It is clear that even assuming the minimum 
value of m found in this analysis, a mild to non-evolving 
merger fraction (m < 1.5) with redshift is ruled out at 
the ^ Act confidence level. 

4.1. Potential Biases 

In many previous merger fraction studies, any evolu- 
tion measured in the merger fraction with redshift typ- 
ically depends strongly on the either the lowest redshift 
{z It 0.2) data point where volume effects come into play, 
or the highest most incomplete redshift bin. Although, 
like many other works our lowest and highest redshift 
bins suffer biases or incompleteness, we still find evolu- 
tion of the merger fraction of at least m ~ 2 over our 
most secure and complete redshift range (0.3 < z < 0.9). 

The interaction fraction for the lowest redshift bin 
(0.1 > z < 0.25) is elevated compared to the GIF found 
between 0.25 < z < 0.6. This is likely related to two fac- 
tors. The first is the larger angular size of low redshift 
galaxies making tidal features more easily detected. The 
second and likely more dominant factor is that the lower 
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Fig. 4. — Mass limited galaxy interaction fraction as a function 
of redshift from the CFHTLS-Deep Survey (circles). The colored 
lines represent the best (1 -|- 2)™ fits when various data points 
are included in the analysis. When all points between 0.1 < < 
1.15 (with stellar masses > 10^'^Mq) are included the best fit of 
(l-t-z)™ is m = 2.25 ±0.24 (blue, solid). When the lowest redshift 
point is not included, m = 2.8±0.25 (red, dash-dot), only the high- 
z point removed, m = 1.95 ± 0.2 (green, dash-dot-dot), and when 
both the low and high-2 points are not included m = 2.56 ± 0.24 
(purple, long d as h). The b lack dashed line shows the result of 
ILin et al.l l(2003 'l: ILotz et al. l ( 2008a) (m = 0.51 ± 0.28) who find 
mild to no evolution in the close pair fraction. The error bars are 
derived using Poisson statistics, while the horizontal errors come 
from the uncertainty in the photometric redshift. The vertical 
dotted line represents the 75% completeness limit. The open circle 
notes t he d ata point suffering from the selection bias discussed in 
Section im 

redshift bins are more susceptible to cosmic variance as 
shown in Figure [5] and discussed further in Section 14.21 

Another potential bias to consider is that as we probe 
higher redshifts, z > 0.9, we begin to more closely ob- 
serve the UV (u*-band) which is dominated by massive 
young O and B stars. Therefore, star formation in tidal 
tails and bridges may be more easily visible resulting in a 
higher GIF at larger redshifts. Since the M*-band images 
are not a comparable depth to the I'-band data, a simu- 
lation similar to that outlined in Section [3.41 was run on 
the sub-sample of lower redshift tidal tailed galaxies but 
with a 1 mag enhancement in surface brightness. The 
resultant GIF increased by ^ 1% at z = 1. It is therefore 
unlikely that this potential bias is solely or even largely 
responsible for the higher GIF at z 1, especially con- 
sidering that the completeness test estimates that we are 
missing a much larger fraction, 15%-30% of tidal tailed 
galaxies, at z ~ 1 due to cosmological effects. Despite 
these potential biases, our analysis suggests that the evo- 
lution in the GIF is a real physical effect, meaning their 
were more galaxies undergoing mergers at earlier times 
out to z ~ 1. 

4.2. GIF: Cosmic Variance 

A benefit of the CFHTLS-Deep survey is that it is 
spread over multiple fields allowing for the study of how 
cosmic variance can affect measurements. In Figure [U 
the GIF is presented separately for the Dl and D2 fields. 
Both fields show an evolving GIF with redshift, how- 
ever the steepness of the evolution varies, stressing the 
importance of studying multiple large fields. The jack- 
knife errors, representing the impact of cosmic variance. 
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were calculated by considering the difference in the GIF 
between the Dl and D2 fields (see Figure [5] inset). At 
lower redshifts, where the volume being probed is less, 
we find the GIF is affected more strongly by cosmic vari- 
ance, while at z > 0.6 little variation in the GIF is seen 
between the fields. 
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Fig. 5. — Galaxy interaction fraction as a function of redshift 
from the CFHTLS-Deep Survey for the Dl (filled squares, black) 
and the D2 field (filled circles, red). Cosmic variance dominates 
the errors of the GIF at z < 0.5; however, at higher redshift the 
Poisson errors are larger. Inset panel shows the difference between 
Jackknife (JN) and counting statistic (CS) errors. 



4.3. GIF: Comparison with Previous Works 

A key consideration when attempting to compare re- 
sults from previous merger rate studies is the similarity 
of the objects being probed, i.e., comparable luminosity, 
or stellar mass limits as well as the selection technique. 
Numerical results of the merger fraction, using morphol- 
ogy should not be compared directly to a close galaxy 
pair fraction, since merger selection techniques have very 
different observability time-scales. Direct comparisons, 
however, of the galaxy merger rate derived with different 
selection techniques, can be made since the close pair 
or merger fractions are normalized by the appropriate 
timescale. In Section 15. 1[ we compare our merger rate 
with those of previous studies. It must be noted that 
the merger timescale can vary widely depending on the 
initial merger parameters and remains a large source of 
uncertainty. Therefore, a comparison of merger fractions 
is warranted, since they are proportional to the evolution 
of the merger rate. 

In this section, we compare our results of the merger 
fraction and its evolution with redshift with other mor- 
phologically based works. Since we use a morpholog- 
ical approach in identifying interactions a first-order 
comparison would be studies that also consider the vi- 
sual appearance of a galaxy whether it be qualit ative 
or quantitative, like t hose of ILe Fevre et all (2000) and 
iGonsehce et al.l (|2003| ) to name a couple. Figure [6] com- 
pares our work with other morphologically based esti- 
mates of the merger fraction. The only other study which 
limited their sample by stellar mass, rather than lum i- 
nosity in this redshift range is IGonsehce et"all (f2003h . 
At intermediate redshifts we find a merger fraction of 



6%± 0.4%, in good agreement with IGonsehce et al.l 
(|2003f ) who find 7% ±5%, and consistent results at z ~ 1. 
Reasonable comparisons can be made with other studies 
since we are also probing similar luminosity ranges (see 
Section 23]). A major obstacle that has plagued merger 
rate studies is small sample sizes inspecting only a few 
tens to a hundred galaxies, resulting in large (5% — 15%) 
uncertainties making it difficult to confidently measure 
both the merger fraction and its evolution with redshift. 
In this study, we inspected more than 25,000 galaxies, a 
factor of > 10 more than previous morphological investi- 
gations. We find that within the uncertainty most of the 
studies shown in Figure [B] are consistent with our GIF 
measurements. 
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Fig. 6. — Plot of the mass limited GIF as a function of red- 
shift from the CFHTLS-Deep Survey using the primary inter- 
acting sample (filled circles, black) and the secondary sample 
(squares, red) that include the primary merger sample and galax- 
ies with "short" tidal tails, and double nuclei. Blue circles high- 
light the interaction fraction corrected for complete ness. Ob- 
servatio ns of the me rger fraction are co mplied from ILotz et al.l 
(!2008al, blue stars), IBridge et alj |20071, green, inverted trian- 
gles), De Propris et al. (2dD 7|, bl a ck x), jCass ata e t a l. (2005, pur- 
ple diamonds), Consclicc et al.l I I2003I . orange, triangles), and 
Le Fcvrc et al. (2000, cyan, open squares). The black lines outline 
the best fit of form (l+z)™ for the various samples, m = 2.25ib0.24 
for the primary merger sample (dashed), while the GIF corrected 
for completeness (dotted line) has a stepper fit of m = 3.31 ± 0.18. 
The best fit for the primary sample with the inclusion of proba- 
ble interactions is m = 2.14 it 0.17 (dash-dotted). All fits include 
the lowest redshift points of this stu dy ( open circle/square) which 
suffer from a selection bias (Section 14.111 . strengthening m = 2.25 
as a lower limit. It is clear that a non-evolving merger fraction is 
inconsistent with the CFHTLS-Deep data. 



This work has presented a statistically secure lower 
limit of the interaction fraction for galaxies with stel- 
lar masses M■^, > 10^'^Mq as a function of redshift. The 
fraction of galaxies undergoing tidal interactions is surely 
higher than the rate reported here, as this study is lim- 
ited to higher surface brightness features, and only in- 
cluded the most confident tidal features in its analysis. 

In Figureini we recalculated the GIF applying the com- 
pleteness correction estimated in Section lX^ Recall that 
this correction factor accounts for the number of mergers 
not identified as a function of redshift due to cosmological 
effects. Evolution of the galaxy merger fraction is strik- 
ingly apparent. We find that the best fit of (GIFo,m) 
is (0.015 ± 0.002,3.31 ± 0.22). Also shown in Figure [B] 
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is the GIF with the inclusion of less confident signatures 
of tidal interactions or mergers, such as galaxies with 
"short" tidal tails (tail length less than diameter of host), 
and galaxies with a double nuclei and no secondary signs 
of an interaction. This resulted in an additional 767 in- 
teracting galaxies being included in the analysis. The 
resultant GIF is on average ^ 5% higher than our lower 
limit, and its evolution with redshift (m = 2.14 ± 0.17, 
compared to m = 2.25 ± 0.24) traces the lower limit GIF 
very well over the full redshift range. Although these 
classifications were deemed to be less confident, they are 
likely dominated by galaxies which are truly merging as 
there is reasonable agreement with the completeness cor- 
rected GIF at high redshift. 

Overall we find that between 0.1 < z < 1.2, our value 
of m {m = 2.25) is in good agre ement with m = 2.2 ±0.3 
found by [C assata et al.' (2005*), m = 2.1 2 repo rted by 
[Bridge et al.. (,2007.) and Conselice et all (pOOl . The 
merger fraction acco unting for incompleten ess is in bet- 
ter agreement with iLe Fevre et al.l ()2000[ ) which find 
m = 3.4 ± 0-6. Our results do differ, however, from 
ILotz et all (|2008aD (m < 1) which utihzed the DEEP2 
survey. There are some factors that may explain the 
discrepancy aside from the choice of luminosity versus 
stellar mass. First, the G — M20 structural parame- 
ters have a significantly short observability timescale (0.2 
Gyr) compared to our visual classification (0.8 Gyr), 
which could account for the lower merger fractions found 
at some redshifts. Second, mergers were identified us- 
ing the Gini — M20 region calibrated by the location o f 
low redshift ULIRGs on this diagram ()Lotz et al.ll200l . 
iKampczvk et al.l (|2007[ ) found that only 6% of visually 
cla ssified merg ers lie in the region of the Gini— M20 used 
bv'L otz et aIT(|2008aD . Lotz et al. 2008 have also found 
through simulations that when a merger has a single nu- 
cleus the Gini— M20 merger classification has a detection 
efficiency of 46%. Third, the mild merger evolution re- 
ported by ILotz et aTl (|2008af ) is strongly tied to the two 
lowest redshift bins which suffer from the largest errors, 
while the merger fractions with the smallest uncertainty 
are in better agreement with our findi ngs. Finally , cosmic 
variance may also have an impact, as lLotz et al.l ()2008al ) 
uses a single field on the sky. 

As discussed earlier, even assuming the minimum value 
of m found in this analysis, a mild to non-evolving merger 
fraction (to < 1.5) is ruled out at the ^4(7 confidence 
level. 

4.4. GIF: Interaction Classes 

The large area of the CFHTLS-Deep Survey allows 
us to further explore the fraction of galaxies interact- 
ing at various stages of the merger process and how 
they evolve with redshift. Figure [7] illustrates the GIF 
for the different merger classifications. The sample was 
again restricted to be brighter than i' <21.9 to en- 
sure confident classifications and to have stellar masses 
logMi,{MQ) > 9.5 to probe similar mass galaxies at low 
and high redshifts. We find that all interacting galaxy 
classes show at least some evidence of evolution with red- 
shift. 

The specific values of the interaction fractions for each 
merger class most likely vary from one another for two 
reasons. First, each classification represents a snap-shot 
of a stage of the merger process, and various signatures 
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Fig. 7. — Galaxy interaction fraction as a function of redshift 
for different interaction classes. Short tidal tails are shown in dia- 
monds (red), medium tails squares (black), long tails circles (blue), 
bridged close pairs upward triangle (orange), galaxies with a double 
nuclei and a tidal tail downward triangle (green), and double nuclei 
with no tail star (violet). The error bars are derived using Poisson 
statistics, while the horizontal errors come from the uncertainty in 
the photometric redshift. The vertical dotted line represents the 
75% completeness limit. A stellar mass limit of 10^-^ (Mq) has 
been imposed. 

have longer durations. For example, a galaxy with a 
tidal bridge can be identified at any point after the first 
encounter, and until the nuclei coalesce, while a galaxy 
with a visible double nuclei and a tidal tail may be visible 
for a shorter timescale as the tidal tail may have faded 
below the SB limit resulting in this classification found 
less frequently. A second factor relates to the resolution 
of the images, as ground-based imaging of very late stage 
mergers, with close double nuclei may be blurred into a 
single nucleus again reducing the number of late stage 
mergers identified. 

4.5. Interaction Fraction: Luminosity Dependence 

After estimating a lower limit for the galaxy interac- 
tion fraction, we now explore the sensitivity of the results 
to various galaxy properties. Many previous merger rate 
studies have fo und the merger rate to depend on opti- 
cal luminosity (iPatton et al.lll997l 120001: fConselice et al.l 
120031: iXu et al.ll2004[ ). In the following sections, we ap- 
ply various optical luminosity limits to our sample and 
re-analyze the GIF to explore any possible luminosity de- 
pendence. Figure [8] shows the absolute g'-band magni- 
tudes (Mq) for the interacting (red) and non-interacting 
(black) populations in the combined Dl and D2 fields. 
It also outlines the Mq limits imposed to aid in a more 
complete comparison of galaxies at high and low redshifts 
with and without luminosity evolution. 

We chose a minimum {Mmin) luminosity limit, Mq < 
—20, which is a few tenths of a magnitude brighter 
than Mq to balance the completeness at high redshift 
with probing Mq as closely as possible. The sample 
was divided into "bright" (Mq < —21.0) and "faint" 
(—21.0 < Mq < —20.0) galaxies to study the impact op- 
tical luminosity may have on the frequency of observed 
galaxy interactions. Figure [H] presents the GIF as a func- 
tion of redshift for the two luminosity ranges. A clear 
dependence of the GIF on AIq is evident out to z ~ 0.7, 
after which bright (blue points in Figure |9]) and faint 
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(red points) galaxies have statistically similar interaction 
fractions. Within the optical luminosities outlined above, 
more luminous galaxies show tidal signatures more fre- 
quently than less luminous galaxies. It must be consid- 
ered, however, that this result may be either entirely or 
in part due to the ease of identifying tidal tails in brighter 
galaxies, rather than a true increase in the frequency of 
interactions in more luminous galaxies. The level of this 
effect is difficult to quantify and requires deeper images 
of these lower luminosity galaxies to see if tidal features 
are evident. 




-ii lftj- . , I , , , I , , , I , , , I , , , I , , , I I 

0.0 0.2 0.4 0.5 0.8 1.0 1.2 200 

Redshift (z) N 

Fig. 8. — Absolute g'-band luminosity as a function of redshift for 
galaxies classified as non-interacting (open circles, black) and those 
undergoing an interaction (filled circles, red). The dashed lines de- 
fine the lower absolute magnitude limits considered in the following 
analysis section; (dashed) no luminosity evolution, (dashed-dotted) 
includes luminosity evolution Q = 1. The histogram on the right 
shows the absolute g' magnitude distribution for the interacting 
galaxy sample. 

Luminous galaxies, below z < 0.6 have interaction frac- 
tions up to 2 times higher than the sample of lower lumi- 
nosity galaxies. Overall, we find that between 0.2 < z < 
0.7 the GIF for more luminous galaxies remains fairly 
constant, while the GIF for fainter galaxies shows mild 
evolution in this redshift range. There is evidence of an 
increasing GIF at z > 0.8, suggesting that interactions 
and merging may play a larger role in galaxy evolution at 
higher redshifts. Although this redshift bin contains '--^60 
mergers, caution should be taken with its interpretation 
as it does suffer from incompleteness, particularly in the 
lower luminosity sample. At this point, we have assumed 
that galaxies in the past and those locally have similar 
optical luminosities. We now proceed with a scenario 
that includes luminosity evolution. 

4.6. Interaction Fraction: Luminosity Evolution 

When a static luminosity limit is considered for a range 
of redshifts, one is essentially assuming that the mass- 
to-light ratio (M/L) of galaxies is the same over that 
redshift range. In this section, we explore the impact that 
luminosity evolution may have on the GIF measurement. 

There is considerable controversy in the literature as 
to how galaxies evolve at z < 1. Although it is agreed 
that evolution does occur ()Lin et al.lll999HCarlberg et al.l 
1200 ID . at what magnitude is under debate. At the very 
least galaxies will evolve passively as their stellar popu- 
lations age, resulting in a gradual fading of their optical 



light. Hence, at higher redshifts when the galaxies are 
younger we would expect them to have higher mean lumi- 
nosit ies. Detailed luminosity function studies ()Lin et al.l 
|1999() showed that a luminosity correction can be ap- 
plied to a galaxy at redshift z, using the expressi on, Qz, 
where Q is typically tak en to be equal to 1 tLin et al.l 
IT99llPatton et al.ll200l) . 

Luminosity evolution has been c onsid ered in some close 
pair studies. iPatton et all (|2002f ) and iLin et al.l ()2004f ) 
found that Nc (the number of companions per galaxy) 
can be significantly affected by the inclusion or exclu- 
sion of luminosity evolution. When luminosity evolution 
is corrected for, fainter galaxies at lower redshifts are in- 
cluded, resulting in higher pair statistics, and in turn less 
evolution of the pair fraction with redshift. 

To explore the effect luminosity evolution may have 
on the interaction fraction we have adopted Q = 1 and 
repeated the GIF analysis described in the preceding sec- 
tion, when Q was assumed to be zero. We find a < 1% 
difference in the GIF with and without the inclusion of 
luminosity evolution. Unlike close galaxy pair's studies 
which typically find an increase in the number of pairs at 
lower redshift when luminosity evolution is applied, we 
find a slight decrease (< 1%) in the GIF , with Q = 1. 
This could merely be caused by the increased difficulty 
in identifying tidal tails in fainter galaxies, or perhaps 
there is an intrinsic dependence of the GIF on optical 
luminosity as implied by the "bright" and "faint" sub- 
samples. 

In contrast to ILin et ahl ()2004[ ) who found a factor of 
4 difference in m with Q = to Q = 1, we find no 
significant deviation in the evolution of the GIF with 
or without a luminosity evolution correction (Figure ^ . 
This discrepancy may be an effect of the different se- 
lection techniques. As you probe further down the lu- 
minosity function, the number of galaxies increases per 
unit volume, which could in turn increase the contamina- 
tion from projection effects and unbound pairs. With the 
tidal tail methodology, probing fainter magnitudes makes 
the identification of mergers more challenging. However, 
the intrinsic close pair fraction or merger fraction may 
also be a function of luminosity. Measuring the pair or 
merger fraction as a function of luminosity is highly prob- 
lematic due to the numerous selection biases and contam- 
ination effects. Disentangling these effects is challenging, 
which is why we have chosen to focus this study on a mass 
selected sample rather than a luminosity one. 

4.7. GIF: Mass Dependence 

Probing the merger fraction as a function of stellar 
mass provides insight into how galaxies build up over 
time. In hierarchical assembly, the leading theory for 
structure formation in the universe, large structures form 
from the merging of smaller structures. This implies that 
the most massive galaxies form latest in the history of 
the universe. Alternative theories postulate that mas- 
sive galaxies can form through rapid collapse, over short 
timescales in the early universe. In this section we aim to 
address the merger histories of low (Af, = 10^ — IO^^Mq) 
and high {M^, > 10^'^-'^ Mq) mass galaxies and the impli- 
cations this has on driving galaxy evolution. 

Using the five-band optical photometry in the 
CFHTLS, we estimated the stellar masses of our sources 
through template fitting (see Section [Z!2l for details) . Fig- 
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Fig. 9. — Galaxy interaction fraction as a function of redshift, as- 
suming no luminosity evolution, Q = (top) and Q = I (bottom). 
The GIF is not strongly dependent on the assumption of luminos- 
ity evolution. At lower redshift {z < 0.6), more luminous galaxies 
(squares, blue) have a GIF up to 2 times that of lower luminos- 
ity galaxies (downward triangle, blue). The error bars are derived 
using Poisson statistics, while the horizontal errors come from the 
uncertainty in the photometric redshift. The vertical dotted line 
represents the 75% completeness limit. 

ure [5] details the stellar mass distribution as a function 
of redshift for our i' < 21.9 sample of classified galax- 
ies. Previous merger rate studies derived from obser- 
vations and theoretical models have suggested that the 
merger fraction not only varies with optical luminosity 
as confirmed in the previous section but also depends 
on stellar mass. It has been suggested that brighter, 
massi ve galaxies have the highest m erger fracti ons at 
z > 1 (jConselice et al.ll2003HConselice 2006; Mallcr_eLalJ 
|2006( ). In order to investigate how the GIF is affected by 
stellar mass, we divided our i' < 21.9 sample into four 
mass ranges; low masses (10^ ° < M^, < 10^-^ and 10^'^ < 
< 10^°), intermediate mass (10^° < M* < 10^° ''), 
and high-mass galaxies {M^, > lO^'''^). The interaction 
fraction is calculated for each stellar mass range. The 
GIF statistics are described in Tableland plotted in Fig- 
ure [101 AH stellar mass ranges imply a similarly evolving 
GIF with redshift. It should be mentioned that the el- 
evated merger fraction at z < 0.6 for massive galaxies 
and z < OA for intermediate mass galaxies (open circles 
in Figure [T0| is likely the same selection bias noted in 



Section O 

A key result of Figure [10] is the indication that lower 
mass galaxies (Af, ^ 10^ — lO^^M©) have a higher merger 
fraction (by ^ 2% — 5%), compared to more massive 
galaxies (M* > IO^^-'^Mq) at z^l. Since previous stud- 
ies at higher redshifts {z > 1.5 — 3) have reported the 
opposite trend, finding the merger fraction of lower mass 
galaxies less than that found in more massive systems 
(Conselice et al. 2009) our result suggests (with moder- 
ate significance) a transition of galaxy assembly between 
z ~ 1 and 1.5 (i.e. galaxy assembly downsizing). An ad- 
ditional implication of our findings is a potential mech- 
anism for interpreting cosmic downsizin g. As mergers 
are known to trigger sta r formation (Mih os fc Hernquisti 
'1996'; ' Barton et a.l.| [2000. and Section [Bj of this work), a 
higher merger fraction in lower mass galaxies would aid 
in the transition of the dominant cites for star formation 
moving to lower mass systems at lower redshifts. 
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Fig. 10. — Galaxy interaction fraction as a function of stellar 
mass and redshift. The data points report the GIF for four stellar 
mass ranges from Log Af=9.0 to 13, and increase in size with stellar 
mass. The open data points likely suffer from the identification bias 
discussed in Section 2] There is similar evolution of the GIF for 
all mass ranges. There is an indication that lower mass galaxies 
at z < 1 interact more frequently than high-mass galaxies. This 
trend is a possible mechanism for downsizing. 

This trend of a higher galaxy merger fraction in lower 
mass galaxies at z < 1 is in contrast with the recent clos e 
pair study in the GOODS fields of lBundv etall (gOOP). 
A possible explanation for the higher pair fractions (of a 
few percent) found in higher mass galaxies is the stronger 
clustering of massive galaxies, which was not considered 
when the contamination due to projection effects was 
estimated. Another factor to consider is the small sample 
size of 1 89 galaxy pairs, compared to our sample of > 
1100 mergers. That being said the trend found in this 
work is still one of only moderate significance, stressing 
the need for even larger samples to better characterize 
the dependence of GIF on host galaxy stellar mass. 

5. GALAXY MAJOR MERGER RATE 

The role mergers play in the formation and evolution 
of galaxies is largely unknown. The rate in which galax- 
ies merge can affect the mass function of galaxies, and 
is likely linked at some level to the decline of the cosmic 
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TABLE 4 

Galaxy Interaction Statistics: Mass Dependence 



Log Stellar Mass (Mq) 


Redshift 






Interaction Fraction [%] 




9.0 < M < 9.5 


0.10-0.30 


1611 


67 


4.2±0.5 




0.30-0.40 


727 


39 


5.4±0.9 


9.5 <A1< 10.0 


0.20-0.40 


1809 


47 


2.6±0.4 




0.40-0.60 


1564 


67 


4.3±0.5 




0.60-0.80 


1216 


90 


7.4±0.8 




0.80-1.00 


180 


24 


13.3±2.7 


10.0 <M < 10.7 


O.iO-0.40 


•li)Vl 


83 


4.1±0.5 




0.40-0.60 


3205 


128 


4.0±0.4 




0.60-0.80 


3107 


207 


6.7±0.5 




0.80-1.00 


1287 


124 


9.6±0.9 


M > 10.7 


0.40-0.60 


1137 


80 


7.0±0.5 




0.60-0.80 


1528 


86 


5.6±0.6 




0.80-1.00 


1060 


102 


9.6±1.0 



Note. — ^Totai is the number of galaxies within the specified stellar mass and 
redshift range. Njnt the is the number of galaxies confidently classified as interacting 
within a given redshift and stellar mass range. The errors for the interaction fraction 
were derived using counting statistics. 



star formation rate. In this section we discuss the galaxy 
major merger rate derived from the CFHTLS-Deep Cat- 
alog of Interacting Galaxies using Equation [2] 



GIF/T„ 



(2) 



where GIF is the galaxy interaction fraction, and T,,ng 
is the assumed merger timescale. This equation provides 
a measure of the number of mergers galaxy"^ Gyr~^. A 
primary source of uncertainty for all merger rate studies 
is the merger timescale. This observability timescale is 
essentially the length of the merger process over which 
a specific technique (e.g. finding close galaxy pairs or 
using morphology) is able to identify the galaxy as a 
merger. Estimates of merger timescales have been de- 
rived using dynamical friction arguments (jPatton et al.l 
120021) as well galaxy sc ale numerical simulations tracing 
stellar or gas mrticles (|Conselicell20"06t iBeh et al.|[2006at 
llono et aLlMil . 

Recent work by ^Lotz et all ()2008al ) , in which a mor- 
phological analysis was performed on a suite of A^- 
body/hydrodynamical equal mass gas- rich mergers that 
have been processed through a radiative transfer code 
have provided reasonable estimates (0.2-0.9 Gyr), for the 
timescales probed by close pair methods and quantitative 
morph ological parameters , such as G — M20 and asym- 
metry. ILotz et all (|2008a|) shows that quantitative mor- 
phological classifications based on G — M20 are sensitive 
during the first encounter and final merger stages for 
gas-rich equal mass mergers, but do not identify many 
interacting galaxies between these two stages, resulting 
in an observability (or merger) timescale range of 0.2-0.6 
Gyr. Tidal tails can remain visible even after the central 
portion of the galaxy exhibits a uniform or symmetric 
appearance, therefore the timescale our technique is able 
to detect mergers is longer than that of G — M2q. 

In order to estimate the merger time-scale probed by 
visually identified mergers based on tidal features, we uti- 
lized iV-bod y simulatio n s of g alaxies undergoing mergers 
descri bed in iConselicd (120061) w hich employs the mod- 
els of IMihos fc HernauistI (|1996h and IMihosI (I200I . We 
carefully considered the duration that a galaxy encounter 
would exhibit the tidal features used in this work to iden- 



tify interacting galaxies. We visually inspected snapshots 
of a simulated merger noting the duration in which the 
galaxies would be classified as "interacting" according to 
our criteria. Ultimately, we estimate the timescale being 
probed by strong visual tidal features to be 0.8 ±0.2 Gyr. 

Using Equation [21 the merger rate for galaxies with 
stellar masses > 10^'^ Mq was derived and shown in Fig- 
ure [TT] as a function of redshift. The large uncertainty in 
the merger timescale reflects the various possible merger 
scenarios (i.e., large mass ratios can extend the merger 
timescale). 

We find an average merger rate of Rmg ^^0.065 mergers 
gal~^ Gyr~^ between 0.1 < z < 0.7, which increases 
to 0.24 mergers gal^^ Gyr^^ at z ~1.0. The merger 
rate evolves with redshift as (1 -I- 2;)2-25±o.23_ -yyhen both 
the Poisson errors in the GIF and the uncertainty in the 
merger timescale are included m = 2.33±0.72. 

Using our GIF results and the merger timescale de- 
rived above, we calculate a lower limit for the interac- 
tion history for typ ical galaxies in our sample. Equation 
[3] ( Consclici I2006D shows that by integrating the frac- 
tion of galaxies undergoing an interaction divided by the 
merger timescale one can obtain the number of interac- 
tions an average galaxy undergoes between two points in 
redshift space. 



r""^ GIFiz) , ,GIFa^, , dz 

n^nt^ I ' dt= tH(^^)(l + z)™-i 



^mg 



E{z) 
(3) 

where tn is a Hubble time, GIF(z) is the galaxy inter- 
action fraction at a given redshift, GIFq is the GIF at 
z -0 , and E{z) = {VIm{1 + zf^iC)''^/'^. A power-law 
increase for the interaction rate was assumed, as it is well 
fit by the data. 

Based on the above equation, and assuming m = 
2.25,GIFo = 0.0215 (the best fit to our data outlined 
in Section |4]) , we find that a galaxy with a stellar mass 
> 1O^'^M0 (average mass ratio > 4 : 1) experiences ■~0.6 
mergers from z = 1.0 to the present day (see Figure [T^ . 
Also shown in Figure [12] is the number of major mergers 
since z ^ \ derived using a range of merger timescales. 
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Fig. 11. — Merger rate as a function of redshift in units of 
mergers galaxy"^ Gyr~^. The filled black squares represent the 
merger rate derived using interacting galaxies with stellar masses 
> 1OS'5{M0). The best fit to our CFHTLS data (solid line) of the 
form (1 + 2)™ finds m = 2.33 ± 0.72,GIF(0)= 0.027 ± 0.003. Com- 
bining the CFHTLS-Deep data with the other works noted below 
results in a best fit with m = 2.83ib0.29 (dashed line) and m = 2.43 
when a 2.0% statistical error is assumed. All fits include both the 
uncertainties of the merger fractions and merger timescale. The 
red points show the observed merger rate derived from close pair 
fractio ns as a function of redshift, from [Patton ct al. (2002, filled 
stars), [Lin et al.' (2004, open upward triangles), Xu ct al. (2004, 
open star). Bell ct al. (200 6H l i ned s tar) , [De Propris ct al. (2007, 
filled square), Kartalte pe et al.l l|2007l . filled upward triangles), and 
[Bundy et al. (2009, open diamonds). The merger rate derived 
from mor phological stu dies as a function of reds hift is shown in 
blue f rom lLe Fevre et al. (2000, op en dia monds), IConselice et al.l 
(|2003, filled circle), fCassata ct jiD (12005!, filled downward trian- 
gles). Wolf et al. (2005, open upward star), Dc Propris et al. (20Q3, 
open star), Scarlata et al. (2007, "x"), JBridge ct al. (2007, filled 
triangle), Lotz ct al. (2008a, lined stars), iConsclicc et al. (20(3, 
open circle), LoDcz-Saniuan ct aT] (|2009l . plus sign), Jogec ct all 
(|2009l , open upward triangle), and this work (filled black squares). 
The assumed merger timescale for merger fractions derived us- 
ing CAS o r concentration, asymmetry was 0.9 it 0.2 Gyr from 
I Lotz et aLl ij^OOSb) and Conseiice ct ahj I I2009I1 . The timescale over 
which mergers selected via G — M20 is assumed to be 0.4 ±0.2 Gyr, 
close galaxy pairs 0.2 ± 0.1 Gyr HLotz ct al. 2008^), and this work 
0.8 ±0.2 Gyr. 

5.1. Merger Rate: Comparison with Previous Works 

As discussed in Section 14.31 a direct comparison be- 
tween close pair fractions and morphological based 
merger fractions is only robust when the measurements 
are converted into a merger rate through the normaliza- 
tion of the appropriate observability timescale. Figure 
[TT] presents the merger rate derived from the CFHTLS- 
Deep survey, and compares it with previous st udies. The 
assume d merger timescales were taken from iLotz et "aLl 
(|2008b[ ). where G — M20 has an observability timescale 
of 0.4 ± 0.2 Gyr, asymmetry (CAS) 0.9 ± 0.2 Gyr, and 
the close pair selection, 0.2 ±0.1 Gyr. 

We find that the merger rate evolves as m = 2.33±0.72, 
GIF(0)= 0.027 ± 0.003. Combining the CFHTLS-Deep 
data with the literature (Figure [TT]) results in a best fit 
with TO = 2.83 ± 0.29, GIF(0)= 0.022 ± 0.003 (dashed 
line). Since the fit is driven by data with small statis- 
tical errors we also performed a fit assuming a standard 
2% error, resulting in to = 2.43. All fits include both 
the uncertainties of the merger fractions and the merger 
timescale. There is striking overall agreement in the liter- 
ature supporting an increasing merger rate from the local 
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Fig. 12. — Interaction history, or the number of interactions an 
average galaxy in our sample has undergone since z ~1 (solid, 
black line) . The coloured lines represent the same measurement but 
assume a different merger timescale, as stated in the plots legend. 
The assumed evolution of the merger rate was (1 + ^)2.24f0.24^ 
The red (dash-dotted) line is presented for comparison purposes 
assuming little evolution of the merger rate. 

universe to z = 1.3. In the previous literature, it is of- 
ten suggested that the various selection techniques (e.g., 
close pair versus morphology) are a source of discrepancy 
in terms of measuring the evolution merger rate evolu- 
tion. Figure [TT] shows that studies in which the close 
pair method was used to identify galaxy mergers (red 
data points) generally agree with morphologically based 
methods (blue points). 

When considering the value of to derived by close 
pair studies o ur fin dings are in g ood agreemeii t with 
IBurkev et afl (IT991 to = 2.5 ± 0. 5), 'Patt on et all (pOOl 
TO = 2.3 ±0.7) and lLe Fevre eFai . (2000 , to = 2.7ib0.6) 
The l argest close pair sample to date iKartaltepe et al.l 
(|2007f) uses the 2 square degree COSMOS field, within 
which our D2 field is embedded. The value of to derived 
using close pairs in COSMOS is mildly steeper than that 
found in our full sample at to = 2.8±0.1, but is con- 
sistent with our findings especially when we exclusively 
use the D2/C0SM0S field (see Figure [S] showing cosmic 
varian ce). It should also be noted that lKartaltepe et alj 
(12001 use photometric redshifts to select galaxies sepa- 
rated by < 20 kpc, resulting in a larger contamination of 
non-merging pairs comparing to spectroscopic samples. 
Since the accuracy of their photometric redshifts evolves 
as 0.03(l + z), close pair fractions at higher redshift suffer 
larger contamination which would drive the evolution of 
the close pair fraction higher, which may explain their 
steeper value of to. 

Our result s do di ffer however from the close pair study 
of iLin et al.] (I2004D (m = 1 . 08 ±0 4) and morphological 
investigat ion o f ILotz et al.l (|2008a| ) (to < 1; as discussed 
in Section lO]) both of which utilized the DEEP2 survey. 



The merg e r rate derived using the close pair fractions of 
ILin et alT ()2004[ ) is consistent with the merger rate val- 
ues found in this work and many others; however, they 
interpret their results as evidence for a fiat merger rate 
evolution. This disparity mainly ster ris from the inclusion 
of low er redshift pair fractions from lPatton et ahj ()2000l 
l200l) . which are higher than other pair fractions in this 
redshift range, when fitting the equation, (1 + z)"^. An- 
other possible source of the discrepancy is that ILin et all 
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(|2004[ ) find that the value of m changes by a factor of 
^ 4 with the assumed luminosity evolution ((3 = 0—1), 
we find no strong dependence of to on Q f Section I4.6p . 
Finally, the small sample sizes are largely unable to sta- 
tistically decipher between mild and moderate merger 
evolution (1 < m < 2); however, the CFHTLS catalog 
having a sample ~15 times larger is not plagued by this 
specific uncertainty. 

6. STAR FORMATION RATES OF INTERACTING 
GALAXIES 

Interactions and collisions can profoundly affect the 
evolution of galaxies, through morphological transforma- 
tions, mass accretion, and perhaps the most brilliant, 
through induced star formation. Obse rvations of inter- 
acting galaxies such as the Antenn ae (lSchweizeHll982f) 
and simulations of major mergers (jMihos fc Hernquistl 
[1996; Cox ct al. 2006) laoth provide evidence that inter- 
actions can trigger violent starbursts. 

As a first step in studying the impact mergers may 
have on a galaxy's luminosity we explored the Mq dis- 
tribution for four primary interacting galaxy types and 
compared them to the average Mq magnitude for the 
non-interacting field population {Ma^-20.66, Figure [T5)) . 
The four interacting types all show a higher average Mq 
of 0.15-0.5 mag, compared to the non-interacting pop- 
ulation (Mg=-20.81 for medium tidal tails, -20.91 for 
close pairs, -20.96 for long tails, and -21.16 for double 
nuclei with tidal tails), suggesting merger triggered star 
formation. 

Taking a step further we used the five-band optical 
photometry of the CFHTLS-Deep survey and derived 
SFRs for each galaxy in our sample (see Section 12. 2p . 
In Figure [141 we show the average SFR for each inter- 
action class and the non-interacting field population as 
a function of redshift. All interaction classes from close 
pairs with tidal bridges (early stage mergers) to galax- 
ies with double nuclei and tidal tails (later stage) ex- 
hibit enhanced SFRs a factor of 1.5-4 times that of non- 
interacting galaxies. The level of star formation enhance- 
ment also grows with the redshift. At higher redshifts, 
one might expect the average SFR to be larger simply 
because more gas is available. We do find that the av- 
erage SFR for field galaxies (which is largely comprised 
of spiral or disk dominated galaxies (> 60%)) increases 
marginally with redshift from ~1.4MQyr~^ at z ~0.2 to 
~3.2M0?/r~^ at z ■~0.75. For interacting galaxies how- 
ever, we see a factor of 3 growth in the average SFR with 
redshift. This increase could be a result of these systems 
being more gas-rich, allowing tidally triggered starbursts 
to be more significant at higher redshifts. 

Figure [14] also shows that galaxies in a close pair linked 
by a tidal tail (an early interaction stage), have SFRs 
similar to field galaxies, while in later stage mergers 
(galaxies with a tidal tail or double nuclei) show the most 
enhancement, especially at z > 0.6. These observations 
agree with iV-body simulations which suggest that an 
initial starburst occurs after the first encounter, but typ- 
ically a larger burst follows at the end of the merger s e- 
quence (jMihos fc Hernauisti[T99llHopkins et al.ll2005D . 

7. SUMMARY 

Using the deep five-band optical photometry from 
2 square degrees of the CFHTLS-Deep survey, we 
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Fig. 13. — Absolute g'-band magnitude (Mq) distribution for 
galaxies classified as interacting and the field. The dotted line 
represents the average Mq for each interacting class. Top: field 
galaxies (interacting galaxies removed) (second, third, fourth from 
top) show galaxies in a paired system with a tidal bridge, and 
galaxies with medium and long length tidal tails and (bottom) are 
galaxies with a double nuclei and tidal features. The dashed line 
defines the average Mq for the non-interacting field population. 
Luminosity evolution of Q = 1 was assumed. 
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Fig. 14. — Average star formation rates for various interaction 
classes, as well as the field population. Inset is the ratio of the av- 
erage SFR of a particular class compared to the field. The field is 
shown by the x's (black), bridged close pairs by triangle (orange), 
long tidal tails, circle (fjlue), intermediate tails, downward trian- 
gle (green), short tails, open diamond (red), double nuclei, square 
(cyan), and double nuclei with tail(s) are shown by stars (purple). 
The errors bars are the standard deviation of the SFRs. 

developed a new classification scheme to identify galaxy 
mergers based on the presence of strong tidal features. 
We visually classified ss 27, 000 galaxies, compiling the 
largest catalog of interacting galaxies (« 1600) in the 
current literature. With this catalog we examine the 
redshift evolution of the galaxy merger fraction for 
galaxies with stellar masses of M > 10^-^ Mq, and the 
impact merging has on star formation. Our results can 
be summarized as follows. 

1. We find the galaxy merger fraction to be 4.3%±0.3% 
at z ^ 3 and 19.0% ± 2.5% at z 1, implying that the 
frequency of galaxy interactions evolves with redshift as 
(1 -|- 2)2-2^^0-2''. This evolution is a lower limit. A mild 
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to non-evolving merger fraction (m < 1.5) with redshift 
is ruled out at the > 4cr confidence level. This merger 
fraction study has the largest sample to date, a factor of 
Z 10 compared to previous work. 

2. The implied average merger rate between 0.2 < z < 
0.6 is 0.075 mergers gal~^ Gyr~^, which increases to 0.24 
mergers gal^^ Gyr^^ at z '--^ 1. Assuming m = 2.25, a 
galaxy with a stellar mass > IO^'^Mq (and average mass 
ratio 4:1) experiences ~ 0.6 mergers from z = 1.0 to the 
present day. 

3. There is a moderately significant dependence of the 
merger fraction on galaxy stellar mass. We find that 
lower mass galaxies with < lO^^M© are more likely 
to be undergoing a merger than more massive systems 
(Af, > IQ-^'^-'^ Mq). This result is consistent with expec- 
tations of galaxy assembly downsizing. 

4. A compilation of this work with previous merger 
rate studies, presents an evolving merger rate, = 
3?o(l + -2)"*, with TO = 2.83±0.29,»o = 0.022±0.003 (er- 
rors include merger timescale uncertainty). Also noted 
is the general overall agreement between close pair and 
morphological selection techniques with respect to the 
galaxy merger rate and its evolution with redshift. 

5. We find that major galaxy mergers have enhanced 
SFRs a factor of 2 higher than non-interacting galax- 
ies. Late stage mergers show the l argest enhancement, 
consistent with TV-body si mulations (jMihos fc Hernquisti 
[1996; Hop kins et aLl[2005l to name a few). 

Overall we have presented a consistent picture of 
galaxy evolution whereby galaxy interactions occurred 
more frequently in the past, typically triggering bursts 
of star formation. These results show that tidal tails and 
bridges can be a powerful and robust tool in quantify- 
ing the interaction and merger rate evolution of galaxies. 
However it must be stressed that the merger frequencies 
presented in this work are lower limits, and the fraction 
of galaxies undergoing a major merger or interaction is 
indeed higher, especially at larger redshifts. 

Another consideration is that much of the star forma- 



tion in the universe is enshrouded by dust. Optical es- 
timates of SFRs can underestimate the effect merging 
has on triggering star formation in dusty systems (C. 
R. Bridge et al. in preparation). It has been shown 
that dusty, luminous IR galaxies (LIRGs) at z > 0.7 
are the dominant population contributing to the comov- 
ing infrared energy density and r epresent ~ 70% of th e 
star formation activity at z 1 (|Le Floc'h et al.ll2005l ). 
So in order to compile a more complete picture of the 
role merging plays in driving the CSFR it is imperative 
that the processes involved in producing stars in these IR 
bright galaxies be understood. We are exploring the role 
merging plays in LIRGs at 0.5 < z < 1.0 in a compan- 
ion paper (C. R. Bridge et al. in preparation). Clearly, 
more work with larger samples is needed to establish the 
merger rate at redshifts > 1. Large survey areas with 
near-IR observations like those that will be provided by 
HST's Wide Field Camera 3 will provide key data sets. 
At low redshift (< 0.1 < z < 0.3) the merger rate is 
also poorly constrained and will benefit from future large 
area, deep, multi-band surveys like that being done with 
the Large Synoptic Survey Telescope ( LSST). 
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